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Abstract We will study the spectral problem related to the Laplace operator in a singu- 
larly perturbed periodic waveguide. The waveguide is a quasi-cylinder with contains periodic 
arrangement of inclusions. On the boundary of the waveguide we consider both Neumann 
and Dirichlet conditions. We will prove that provided the diameter of the inclusion is small 
enough in the spectrum of Laplacian opens spectral gaps, i.e. frequencies that does not propa- 
gate through the waveguide. The existence of the band gaps will verified using the asymptotic 
analysis of elliptic operators. 
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1. Introduction 

The main goal in this paper is to study the spectral properties of the Neumann- 
Laplacian on the periodic singularly perturbed periodic quasi-cylinder r2(e): 

- Au'{x) = X'u'ix), X G 0(e), (1) 
dnu'{x) = 0, X E dn{e). (2) 

This spectral problem should be interpreted as a singular perturbation of the same 
problem in the straight cylinder Q, = uj x R C R^. 

The quasi-cylinder r2(e) is a periodic set which depends on a small parameter 
e > 0. It is obtained from the straight cylinder 17 by a periodic nucleation of small 
voids with diameter of order e (see fig. la and fig lb, respectively). We remark that 
without lost of generality we have reduced the period to one by rescaling. 

It is known (e.g. (28l|) that the spectrum a of the Neumann-Laplacian in the 
straight cylinder implies the continuous spectrum cTc which is the closed positive 
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Figure 1. A cylinder and a quasi-cylinder with periodic family of voids 



Figure 2. The gap between the first and second bands, while a gap between the second and the fourth 
band is covered by the third one 

real semi-axis [R+, and the point spectrum is empty. In particular, for any real 
/.i > the functions, or oscillating waves, 

give rise to a singular Weyl sequence of the problem operator A at the point E IR+ 
(cf. §9.1 in fH). 

The structure of the spectrum a{e) in the perturbed quasi-cylinder ^{e) is 
much more complicated. Owing to the Gel'fand transform [lO| (see Section 2) and 
the Floquet-Bloch theory (see e.g. I21I ) the spectrum cr(e) of the Neumann- 
Laplacian is endowed with the band-gap structure. Namely, the essential spectrum 
cr(e) = cre(e) is a union of closed connected segments T„(e) C [R+: 

a{e) = U T.(6). 

nGN 

We obtain a spectral gap in the spectrum if there exist an interval in the positive 
real semi-axis which is free of the spectrum, see the left part of fig. 2. However, 
when the segments overlap each other (see the right part of fig. 2), no spectral gap 
opens. 

Here we note that the essential spectrum o"e(e) coincides with the continuous 
spectrum provided none of the segments Tn(e) collapses to a single point, which 
becomes an eigenvalue of the operator ^(e) with infinite multiplicity, and belongs to 
the point spectrum while the discrete spectrum is still empty. A negative answer 
to this possibility is still an open question and the authors can only prove the 
following: For any t > there exists et > such that cTp(e) H [0,t] = when 
e G (0, ei]. We do not pay any attention to this incomplete result and therefore 
always speak in this paper about the essential spectrum. 

In the literature the existence of spectral gaps is mainly investigated for periodic 
media which are infinite in all directions while the coefficients of differential oper- 
ators, both in scalar and matrix case, are usually assumed to be contrasting (see 
e.g. [a, 0, S llOl for scalar problems and [1, [l^] for systems). Much less results 
are obtained for periodic waveguides, which are infinite in one direction only. In 
this case spectral gaps ought to be opened by varying the shape of the periodicity 
cell only which is in accordance with the results in the known engineering practise. 



(3) 
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Figure 3. The two-dimensional dummy of the periodicity cell 



There are two useful approaches to fulfil the goal. The first one utilises the 
asymptotic analysis of the spectrum for the associated model problem in the pe- 
riodicity cell, like we do in the present paper. This approach has been realised 
for the Dirichlet-Laplacian in papers [3, Q, H^j and others. We emphasise that 
our results for the Neumann-Laplacian are completely new and require different 
techniques. 

The second approach is based on the application of sharp parameter-dependent 
estimat es p rovided by Korn's inequalities, usually weighted and anisotropic, see 



e.g 



171 . Il8l . l22l | and others. On one hand, this approach does not require a precise 
description of the dependence of the cell shape perturbation on the small geometric 
parameter, but, in contrast to the first approach, one can only prove the existence 
of one or several spectral gaps without any information on the position and width 
of the gaps. 

Combining of the above-mentioned approaches one could get much more elab- 
orated results by using their advantages. However, we cannot examine here the 
whole spectrum and to our knowledge there does not exist any results of this type. 
Here indeed we discuss only the first spectral gap. 

To prove or to disprove the existence of a gap between the segments Ti(e) and 



T2(e) (cf. fig. 2), we employ the method of compound asymptotic expansions [ij] 
to the associated model spectral problems in the periodicity cell 

zu, = w\e^, zu = ujx (0,1) (4) 

of the quasi-cylinder ([5]). The straight cylinder with the small void 9^ (see the 
two-dimensional dummy in fig. 3b) is nothing but a good example of a domain 
with the singularly perturbed boundary dw^ = dwL) dO^. The asymptotic analysis 
of th e eig envalues for the Laplace operator have been developed at a great extend 



in [8|, [Ij, lla, Il6|, I23h25[| and many others. However, the dependence of the model 
problem on the dual variable r/ G [0, 27r) of the Gel'fand transform brings a serious 
complication to both, the formal asymptotic procedures and the justification since, 
instead of a single spectral problem, we have to deal with a family of spectral 
problems dependent on the real parameter rj S [0, 27r) (see Section 3 and Section 4 
for more details). 



2. Preliminaries and notations 

2.1. The quasi-cylinder and the Gel'fand transform 

The purpose of this section is to explain what do we mean by a periodic quasi- 
cylinder. There are several ways to describe it; but here we start with the straight 
cylinder which is the cartesian product of a bounded Lipschitz domain w C and 
the real line: = w x [R. 

Consider now a fixed bounded smooth domain 9 with the coordinate origin ^ = 
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in its interior. We introduce the family of sets Oi by defining 

ei ={x= (2/1,2/2, z) : (2/1,2/2, z - j) G 0° = {x : ^ = e~\x - x°) G 6}, 

where G za = lo x (0,1). Note that the diameter of the set 6i is of order e. 
The quasi-cyhnder ri(e) is a periodic set which depends on a small parameter 

e > and it will be obtained from the straight cylinder Q excluding the sets 9i 
from it (see fig. la and fig lb, respectively). In other words, the quasi-cylinder can 
be written as 

n{e) = n\\J'9i. (5) 

iez 

We call the set vj{e) = tz7 \ 6*^ in (jH) as the periodicity cell of the quasi-cylinder. 
Note that the cylinder can be viewed as a quasi-cylinder with the straight peri- 
odicity cell as in dH). 

We recall here the definition of the Gel'fand transform [lo| : 

v{y, z) ^ V{y, z, r?) = -|= er'''"v{y, z + j), (6) 



2- 



where (2/,z) G r2(e) on the left, rj G [0,27r) and {y-,z) G tn(e) on the 
right. As it is well known this operator is an isometric isomorphism from 
L2(17(e)) onto L2(0, 27r; ^^(^(e))) and an isomorphism between the Sobolev spaces 
H^{n{e)) and 1^0,271; H^{vj{e))) (see e.g. [13,13). Here L'^{0,2tt- L'^{w{e))) con- 



sists of L^(tz7(e))-valued (complex) L^-functions on [0, 27r]. Similarly, the space 
L^(0, 27r; i?^(tz7(e))) contains those H^{uj{e))-valued (complex) L^-functions on 
[0, 27r]. Finally we have denoted by H^{w{e)) the space of Sobolev functions / which 
are quasi-periodic in the z-variable, i.e. such that function {y,z) i— t- e~^^^f{y,z) is 
1-periodic in z. 



2.2. The relationship between waves and Floquet waves in the straight 
cylinder 

Before going into the asymptotic analysis of the spectral problem in the periodic 
quasi-cylinder we will study the relationship of the standard wave solutions and 
Floquet waves for the Neumann-Laplacian in the cylinder 0. They are the solutions 
of the differential equation 

— Axu{x) = Aii(x), X £ Q , (7) 

supplied with the Neumann conditions on the boundary 

dnu{x) = 0, xedn. (8) 

Applying the Gel'fand transform to ([7]) and ([8]) we obtain the model problem in 
the periodicity cell zu with the dual Gel'fand parameter r] G [0,27r), namely, the 
differential equation 



- A^U{x,r]) = A{r])U{x,r]), {y,z)£zu, 



(9) 
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the Neumann condition on the lateral boundary of 7 = dvj 

dnU{x;r,) = 0, xej, (10) 

and the quasi-periodicity conditions on the opposite ends of the cylindrical cell 
■cu = u X {0, 1): 

U{y, 0) = e-'miy, 1), d,U{y, 0) = e-'^^d^Uiy, 1), y € u. (11) 

Let Ap be an eigenvalue and let Up be the corresponding eigenfunction of the 
model spectral problem ([8|)- (fTT]) . then the Floquet wave can be written as 

n3{y,z)^e''^^'%{y,z-[z]) (12) 

where [z] G Z stands for the entire part of the real number z. Notice that function 
(fT2]) is smooth due to the quasi-periodicity conditions pT]) and it becomes a solution 
of the homogeneous problem ^ with the spectral parameter A = Ap. 
Any function of the type (jl2p must be obtained from the standard waves 

(y,z)^e<^y,(y) (13) 

in the cylinder $7, where C S K is the dual Fourier variable and Vq is an eigenfunction 
of the model problem 

- AyV{y) = MV{y), yeu, d^Viy) = 0,y € dco (14) 

on the cross-section co. Notice that the wave solution (fT3|) with q = 1 and Vi{y) = 
const corresponds to the eigenvalue Mi = and coincides with the wave solution 
([3]) where 

f^ = e 

stands for the parameter in the Helmholtz operator A^. -|- /u. The main difference 
between formulae (I12p and (llSp follows from the role played by the parameters 
77 E [0, 2tt) and C S R. 
We write the wave solution (I13p as follows: 

e<2"l/g(i/) = e*'^^~^™)We**^''~^'^'"^*^^~W''e^'^™'^yg(?/) (15) 

where m = [(27r)~^C] is the maximal integer such that 27rm < (. Now we denote the 
first multiplier on the right of (I15p by e*^'^' and the remaining part by Up{y, z — [z]) 
so that 

Up{y,z) = e'^'e^^^'^Vg{y), (16) 
r] = C- 2TTm e [0,27r). 

Since e^'^™^ is 1-periodic in z, we see that the function Up satisfies the quasi- 
periodicity conditions (fTT]) . In other words, (fT5]) implies the Floquet representation 
(fT2]) of the wave solution with the ingredients p^ . 
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Figure 4. The lowest dispersion curve and the lattice frame made from it 

Relationship (fT6]l helps us to list all eigenpairs {Ap, Up} of the model problem in 
the straight cell w in terms of the eigenvalues 

= Ml < M2 < Mg < . . . (17) 

and eigenfunctions {Vi, V2, V3, . . . } of the model problem ([H]). Namely, there holds 
the relations 

Ag,m{v) = Mg + {rj + 2iTmf , q € M, m € Z. (18) 

We emphasise that the eigenvalues Ag^miv) ^^ow enumerated with two indexes 
so that it is necessary to reorder them to get the monotone increasing unbounded 
sequence 

Ai(??) < A2(r/) < ••• < Afc(77) < ••• . (19) 

The lowest dispersion curve (fTSj) in fig. 4a, corresponding to the smallest eigen- 
value Ml = of the Neumann problem in u, divides into an infinite family of 
smooth finite curves which by above formulae are joined to the lattice in fig 4b. 
In the sequel we are mainly interested in the couple of lowest curves in the lattice 
which are redrawn in fig 5. a, 

A" (,?) := Ao,o(r/) = ^^ A" (r?) := Ao,-i(r?) = (r? - 2vr)2. (20) 

We shall use the notation (j2U|) throughout the paper. 
The variational formulation of problem ([9l)- (llip reads as 

(V,.[/, V^V)^ = A{r]){U, V)^,Ve H^iw) . (21) 

Since for any real parameter rj the bilinear form on the left-hand side of (j2ip is 
symmetric, closed and positive, we can associate with it a self-adjoint positive 
operator in the Lebesque space L^(ro) pj, Thm. 10.1.2]. Moreover, the compact 
embedding H^{w) C Lp'{w) and Theorems 10.1.5 and 10.2.2 in [1| ensure that 
indeed the spectrum of this operator is discrete and forms an unbounded non- 
negative sequence (119p which is constructed from (jl7p . 

2.3. Perturbation of the dispersion curves 

In the perturbed periodicity cell ro^ the both dispersion curves (I20p experience a 
perturbation and, in principle, there can occur two different situations pictured in 
fig. 5.b. and 5.c. First, the dispersion curves [0, 2ii) B t] ^ ^±(^) may still intersect 
so that the spectral segments Ti (e) and T2(e) touch each other and the gap between 
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Figure 5. Unperturbed (a), intersecting (b) and disjoint (c) curves 




Figure 6. Overlapping bands. The dotted lines correspond to the perturbation of the lowest curves in fig. 
5.a. 

them is closed. The second possibihty is that the graphs of the functions 

do not intersect and the area overshadowed in fig. 5c projects into a gap between 
segments Ti(e) and T2(e). In the sequel we will describe when the gap appears 
and indicate asymptotically its position and width. 

Before proceeding with the problem statement, the formal asymptotic analysis 
and the justification, we note that the gap opens only under the condition 

Mi+7r2<M2. (22) 

Otherwise, the segment T2(e), which corresponds to the perturbed eigenvalues in 
(jlSp with q = 2 and m = 0, —1 covers the gap opened between the segments Ti(e) 
and T3(e) (see fig. 6) which are now generated by eigenvalues in (llSp with q = 1 
and m = 0,-1. 

However, the rescaling performed in the very beginning rescales also eigenvalues. 
To see this, assume that the period / of the void nucleation is bigger than one. 
Now the rescaling diminishes the size of the cross-section, i.e. the change of the 
variables {y,z) — )■ {y' , z') = {l~^y,l~^z) yields the shrunken cross-section 

oj' = {y' e : ly' G uj}. 

Simultaneously, eigenvalues of the model problem in a; changes into the eigenvalues 
M[ = and M2 = PM2 of the model problem in the shrunken cross-section to'. By 
choosing / large enough we fulfil the condition (I22p for M[ . In other words, posing 
the small voids 9^ with a sufficiently large period may open a gap in the essential 
spectrum, where as if the period of the voids is sufficiently small first triple of the 
spectral bands surely overlap. 

Assuming condition (j22p we study the spectral problem ([I])- ([2]). Applying the 
Gel'fand transform ^ to problem ([II)-©, we obtain the model spectral problem 
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in the periodicity cell vj^ depending on r/ G [0, 27r): 

- /^W{y, z) = A'{v)U'{y, z), {y, z) e , (23) 
dnU'{x) = 0, X G 7 U de{e) , (24) 
Wiy,0)=e-''^U%y,l), d,U%y,0) = e-'^d,U%y,l), y G ^q. (25) 

As above, using Theorem 10.1.2 [1] we associate with these spectral problems a 
family of unbounded positive self-adjoint operators T'^irj) on L'^{vu^) with discrete 
spectra, i.e., the eigenvalue sequences 

< K\{ii) < Alirj) < ■ ■ ■< AUrj) < • • • ^ oo. 

The immediate objective becomes to construct asymptotics for these eigenvalues 
as e — > +0. 



2.4. Mixed boundary value problem 

Although all the calculations are presented for the Neumann conditions on the 
boundary of the quasi-cylinder 0(e), we also discuss slight modifications, which 
are needed to adapt our calculations to the mixed boundary value problem 

- An'(x) = X^u^{x), X G n{e) , 
dnu'ix) = 0, X G dn{e) \dn = Uj^zdOi , (26) 
u'{x) = o,xedn. 



We emphasise that the asymptotic procedure of [IJ] employed here is not sensitive 
to the type of boundary conditions on the lateral part of the boundary dco x (0, 1) 
of the periodicity cell is employed. However, the change of the boundary condition 
from Neumann to Dirichlet on the surfaces of small voids crucially simplifies the 
asymptotic analysi£]. In this way our paper becomes much more technical than 
(20I I and [2[, where the Dirichlet problem was considered in the regularly and sin- 
gularly perturbed periodicity cells. The particular conclusions on the geometrical 
characteristics of the gap are also different. 



3. The formal asymptotic procedure 

In this section we present a calculation which furnishes an asymptotic formula for 
the first two eigenvalues (j20p . In order to simplify the presentation we start with 
considering the case when the eigenvalues for the model problem are simple, e.g. 
when the Gel'fand dual variable rj is not equal to vr. Notice that the dispersion 
curves in fig. 5a intersect each other just at r/ = vr, and we examine this situation 
in the end of this section. 

In Sections 3.3 and 3.4 we give detailed calculations of the boundary layers and 
a regular terms which lead to a direct formula for the main correction term in 
the asymptotics of the eigenvalues. With this asymptotic formula we are then able 
to consider also the case 77 w tt although at rj = tt where the model problem 
has a multiple eigenvalue. As a result, in Section 3.5 we derive a system of two 



^This does not hold in the two-dimensional case, where the existence of a gap is still fully an open question. 
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linear algebraic equations from which we can find out the perturbation terms in 
the asymptotics of the eigenvalues when ry approximately equal to it. 

To make the necessary conclusion for the case r/ ~ vr, we will apply the trick 



proposed in 20(] and introduce the additional deviation parameter -0 in the repre- 
sentation 

= vr + • (27) 

With the help of this new variable '0 G [— ■i/'OjV'o] we can write the asymptotic 
formulae at the ce^-neighbourhood of the point rj = it. 

We emphasise here the following difference in the notation: In Sections (|3.ip - (l3.4p 
we set 

A%{v) = Alirj), At{r^) = A^(r/), r? G [0,7r), (28) 
A;(r/) = A^(r7), Ai(ry) = AUrj), rj G (7r,27r), 



while in Section (13. 5p 

AXirj) = A^rj), Ai(r/) = Al{r^), t? = n. (29) 

The evident difference in the notations (|28p and (1291) helps to keep similar asymp- 
totic formulae in both cases. 

Calculations in Sections 13.1113.41 are given for the Neumann problem ([l])-([2]) in 
such a way that they can easily adopted for the mixed boundary value problem 
(i26l) . Namely, the eigenpair of the problem (fT4|l {Mi = 0, Vi = const} must be 
replaced by the first eigenpair {Mi, Vi} of the Dirichlet problem 

- AyViy) = MV{y), y£io, F(y) = 0, y £ du . (30) 

Notice that the eigenvalue Mi is positive and simple while the corresponding eigen- 
function Vi can be chosen positive inside the cross-section w. 



3.1. Asymptotic expansions 

To describe the asymptotic behaviour of two first eigenvalues Af^ (ry) , we accept the 
following ansatz: 

A^(r/) = AO (r/) + e^A'^iv) + A^(,?), (31) 

where Aj.{rj) is the eigenvalue (j20p of problem ([7])- ([9]), A^(?y) is a correction term 
and Af^(ry) a small remainder to be evaluated and estimated. 

According to general asymptotic method in domains with singular perturbations 
of the boundaries (see e.g. [IJ]) the corresponding asymptotic ansatz for eigenfunc- 
tions reads as follows 

UUx;v) = Ul{x;r,) + ex{x)Wl{e-\x - x%v) 

+ e\{x)wUe-Hx - x'');v) + e'u:^{x;r,) + UUx;v)- (32) 

The functions are of boundary layer type and x is a smooth cut-off function 
which equals one in a neighbourhood of the point x^ and vanishes in the vicinity 



October 28, 2011 



2;14 Applicable Analysis BANARU gAPA 



10 

of the boundary dvj. The functions are the first two eigenfunctions of the 
unperturbed problem (cf. formula (|16p ). namely 

C/0(x;7?) = e^''^yi(y), Ul[x-r^) = e'^^~^^^'Vi{y), (33) 

where Vi is the first eigenfunction of the Laplace operator on the cross-section, 
normalised in ^^(tu). Clearly, both the functions ([33]) meets the quasi-periodicity 
conditions (jlip . Note that in the case of the Neumann boundary conditions on 50 
the function Vi is just a constant. One readily verifies that since Mi = is the 
first eigenvalue of the Laplacian on the cross-section in the Neumann case, the first 
couple of eigenvalues of the unperturbed problem take the form (j20|) . 



3.2. The boundary layers 

The boundary layer VFj. depends on the "fast" variables ("stretched" coordinates) 

C = (6,6,6) = e-'(^-x°) 

and is needed to compensate for that the main regular term in the expansion 
does not fulfil the Neumann type boundary condition 

dn{^x)u{x) = 0, X G de^. 

By the Taylor formula applied to the function in the "slow" variables x, we 
have 

5„(5)C/l(x;r/) = n(6 • V,[/^(xO) + 0(e), x G dO,. 

Thus Wj. ought to be a solution of the exterior Neumann problem 

A5VFi(6 = 0, (34) 
dni^^Wlii) := GliO = -n(6 • V,C/l(xO;7?), ^ G d9. 

Since Jqq nk{C) ds{^) = for A; = 1, 2, 3, the function G\. satisfies the orthogonality 
condition 



Jae 



Thus, the problem (I34p admits a solution which decays faster at infinity than the 
fundamental solution ^{C) = (4vr|,^|)~^ of the Laplace operator in R^. Moreover, 
the solution has the asymptotic behaviour 



wUO = E Q,k^(0^i^';v) + om~') (35) 



where Q = {Qjk)j,k=i,2,3 is the matrix of the virtual mass for the set 6 C IR'^, which 
is symmetric and negatively definite for any subset 6 with a positive volume (see 
[H, Appendix G]). 
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Remark 3.1. As shown in [26', Appendix G] (see also the representation 

Q = -I3meas3(e) + Q (36) 

is vaUd, where I3 is the 3x3 identity matrix and Q is non-positive. The latter 
is still negative definite in the case meas3(0) / 0; but degenerates for a straight 
crack. 

In order to find the appropriate boundary value problem for the second boundary 
layer term W"^ we take into account a quadratic polynomial in the Taylor formula 
for Ul: 



3 orrO 1 3 ^2rr0 

daUl{x-r^)=Y,a,^{x';r^) + - J] da{x,x,)—^{x';r^) + 0{\x - x^)- 
3=1 ^ j,k=i ^ 

In this way, we conclude that the second boundary layer term should satisfy 
the following exterior Neumann problem: 



The function admits the representation 

WUC) = c±p-' + 0{p-^), p = \C\^oc. (37) 

To calculate the coefficient c± in (|37|) we use the Green formula in the domain 
H/j = Bji{0) \ 9 where B^(0) is a ball of radius R centred at = to obtain 



Ir op Jq0 



dWl 



n dp 
3 Q2rr0 



(0 ds^ (38) 



Notice that on the sphere Sr = dBji{0) we have 9„(g) = dp. 
Applying the asymptotic expansion (f37|) . we see that in ([38 



lim / ^^{Ods.= lim / ( + 0{p-^)] ds. = -Attc. 



'± ■ 



Using the Green formula again but now inside 0, note that the normal becomes 
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inward then, we derive the relation 

j,k=l ■' j=l J 

= AJ7^(x°)meas3(e) = -A^(r?)C/^(x°)meas3(^) . 

In the last equality we also have used the fact that is the eigenfunction of the 
unperturbed problem corresponding to the eigenvalue A^(r/) ( see Q). 
Collecting the above results we finally obtain the formula 

47rc± = -A^(r/)C/^(x°)meas3e. (39) 
3.3. The regular correction 

For the calculation of the main correction term (r/) in the asymptotic expansion 
of A^(?7) we substitute the asymptotic ansatze ()3ip and (I32p into the spectral 
problem (I23|)-(IMI)- By ^ and ^ the boundary layer terms Wl{e~'^{x - 
x^);r]) and W^{e~^{x — x^);r]) have the following behaviour: 

Wl (kx-x%ri]=e' (40) 
+ 0{e^\x -x°\-^), 

Wl{e-\x -xy,rj) = ec±\x - x^\-^ + 0{e^\x - x^\'^) . (41) 

Hence we notice that in (j32p the contribution of the boundary layer terms written 
in the "slow" variables looks as follows: 

e^VF^(e-Hx-x°);r/) + e2w|(e-^(x-j;°);r?) = e^w±{x;7]) + 0{e^\x - x'^]-^) . (42) 

In other words, the main regular correction is of order e^. We write the term ■w±{x) 
as a sum 

w±{x) = w]^{x) + wl (x) , (43) 

where 

^i(-) = - E Q.-^f^(^°) .?"l3 ^ = (44) 

^-^ oxj 4-k\x — x^r \x — x^\ 

j,k=l J I I II 

Now substituting the asymptotic ansatz ([32]) for U^. in the equation (j23|) we get, 
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using that is the solution of unperturbed problem, 

+0{e^\x-x^\-^). 

On the other hand, since A^(r/) is assumed to be an eigenvalue, we also have the 
relationship 

-A,Ul (x; r/) = (A" (r,) + A'^ (7?)) (C/O {x;r,) + {x{x)w± (x) + C/^ (x; r/))) 
+0(e^|x-x°|). 

Thus, setting the coefficients on equal, we obtain the following equation: 

- (A,. + A" (,?)) U'^{x;v) = A'^mUx;v) (45) 
+ {A^ + Al{7])){x{x)w±{x;7])), x€w 

with the boundary and quasi-periodicity conditions 

a„[/4(x;r/) = 0, XG7, (46) 
i74(y,0;r/) = e-*^C/4(y, 1; r/), y G 70, (47) 
d,U:^{y,0;v) = e-'^d.U:^{y,l;v), V e To- 

By the Predholm alternative the boundary value problem ()45p -(j47 p has a solution 
Uj^ if and only if the right hand side of the equation (I45p is orthogonal to the 
eigenfunction in L'^{w). The normalisation ||C/^; L^(g7)|| = 1 then yields the 
expression for the correction term 



A'^iv) = A'±{v) Ul{x;rj)Ul{x;T])dx (48) 



Ulix; 7]) (A^ + A^(ry)) ix{x)w±ix; v,)) dx . 

Recalling that the functions in (j44p are harmonic, we have 

(A, + Al{r,)){x{x)w±{x-r^)) = 0(|x - x'\-^). 

Hence the last integral in (j48p converges absolutely, and therefore it can be calcu- 
lated as a limit of an integral over the set w \ Br{x^) as r — >■ +0. Since the cut-off 
function x vanishes at the boundary of w we obtain by the Green formula that 

Ai(„) . - / (I^a„,.,»,(.; „) - (4.) 

T-->+0 Js V / 

Substituting w± in ([39]) we get 

A'^iT]) = Ki + K2 + Ks + 
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where 



^0 = - 1™ / Ul{x-'q)dnMW^j^{x;ri)dsx, j = 1,2, 
Kj+2 = lim / wi{x; ini)dn(x)U±{x\ rfjds^, j = 1, 2. 



The term Ki is computed by a direct integration: 

j,k=l -^'j '^^+^^8.(0) 



^1 = .11% /3^(„) + ^)^"(^) 4^ 

E Q,.^(.-.)^lirn^ / f/^(xO + .;.)(^) d.. 



^(V.C/°(:EO;r/))TQV,C/0(xO;r/) 



Since dn(x)w\{x\r\) = — 9p (c±p ^) + 0(/9 ^) on Sr(a;'^) (see ([37|) we have 
= lim / W^)^ f 



For the computation of Kj, we write 



w 



^ /5r/0 _ 



.■^ 5xj ' 47r|x — x' 

This yields us 



3 ^fyO /■ ^2 ^jjO 



Sp(0) 

= iv,C/°(x°;r?)TQV,C/0(xO;r/). 
Finally, it is easy to see that 

^^4 = - lim / - — ^—^dn(x)U\(x\r\)ds = ^. 

Collecting the previous results we finally obtain the formula for the main asymp- 
totic correction term 

A'±(??)=F±(r?) 
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where 

F±{v) = V,C/0(xO;77)Tgv,,C/0(xO;r/) 

+ (Ml + (?7 - vr ± 7r)2) |C/^(x°; 77)^ meas3(6') . 

The remainder A^(r/) in the asymptotic expansion (j3ip will be estimated in §4. 
Note that -F+(vr) = -F-(vr) and denote it by Fq, so that 

Fo = V^UI{x°;tt)~^QVJJK^^+ (Mi + vr^) |C/0 7r)|^ meas3(^). (50) 
3.4. Perturbation: 77 ~ tt 

In this section we will study the behaviour of the eigenvalues (rj) when rj is close 
to vr. For the unperturbed problem, i.e. when e = 0, the eigenvalue Aj.(7r) = Mi+vr^ 
has multiplicity two. Due to this change of the multiplicity our analysis becomes 
slightly different. However, we may use the previous calculations. 

In the vicinity of = vr we apply the formula ()27p . In this case the asymptotic 
ansatze (I3ip and (1321) are still valid but their entries depend on the deviation 
parameter ip. Moreover the main regular part of the expansion becomes the linear 
combination of the eigenfunctions and : 

U^(x; vr, V) = af{ij)Ulix; vr) + afitl^)U^ {x; vr) 
= (afWe--+a±(V)e--)yi(y) 

where the coefficient vector a='=(^) = {of (ip) , (ip))^ is to be determined. Without 
loss of generality we may assume that |a^| = 1. By the same argument as in the 
previous section, we have the boundary layers W^and as in (j40p and (|4ip . 
However, in this case they are the linear combinations 

W^'(e;vr,V') = af{ij)W^{^;7T) + af{i;)wUC;7T), j = 1,2. (51) 

Then equation (j45p takes the form 

(-A, - AO(vr))U'±(:E;vr,V^) = A'(^,V')Ui(x;vr,V^) 

- (-A^ - (Ml + TT^)) {xix)w{x; vr, V)) , x£w\ x°, 

where the w(x; vr, ip) is calculated for the boundary layer terms (I5ip according to 
the formulae (1421) and (I43p . The boundary and quasi-periodicity conditions (I46p 
and ()47p turn into following: 

dnU'±{x;7r,i;)=0, x G 7, (52) 
U'±(y, 0; TT, V-) = e-^-U'±(y, 1; vr, V') - ii^e-'^UHy, 1; vr, V'), y S 70 (53) 
9,U'± (y , 0; vr, V') = e'^^a^U'^ (y, 1; vr, ^) - #e"'"9,Ui (y, 1; vr, V') , y G 7o • 

Note that the extra terms on the right-hand side of ()53p is due to the original 
quasi-periodicity conditions (f25]) and the Taylor formula 
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To find A'(7r, V') we multiply equation (f52]l by U^{x;it) and integrate over zu 



?70(x;^)(-A, - A0(7r))U'±(x;7r,V)dx = A'(7r,^)af (V-) (54) 



Ulix; tt) (-A^ - (Ml + TT^)) ix{x)w{x; tt, ^j)) dx. 

The integral on the right can be calculated in a similar way as in the previous 
section and it will provide the formula 



Ulix-^) - (Ml + ^2)) (x(x)w(x;7r,V)) dx = a±(V')Fo + 

+ af (V) (v,[/°(xO;7r)^QV,[/0(xO;7r) - (Mi + 7r2)C/0 (x"; 7r)f/0 (xO; ^)) 



= af(V^)Fo + a±(V)Fi, 
where Fq is given in ()50p and 



i^i = V,?7°(xO;^)'QV,C/0(xO;7r) (55) 



+ (Ml +^2)C/!^(x'^;7r)[/_^(xO;7r)meas3(0). 

The integral on the left hand side of (|54p is calculated by means of the Green 
formula 



Ul [x- vr) (- A, - A°(7r))U'± (x; vr, V-) 



9c 



(U'±(x; ^, V)5nt/0 (x; ^) - a„U'±(x; ^, ^)C/0 (x; tt)) ds,. =: I . 



For calculating this surface integral we take into account the Neumann conditions 
(i52l) for U'_|_ and pUj) for C/^ on 7, the quasi-periodicity conditions (fTTj) and 
to obtain 



- (y, 0; 7r)U'±(y, 0; vr, V) + t/^(?/, 0; 7r)5,U'±(y, 0; vr, V')) dy 
+ / (y, 1; 7r)U'±(y, 1; vr, V') - [7° (y, 0; 7r)9,U'±(2/, 1; tt, ^)) 

= 27rV'af(V) / \Vi{y)\''dy = 2^^4{^). 



7 



In last equality we have used the normalisation condition for Vi and the trivial 
equality \\Ul-L\w)\\ = \\Vr,L^m. 
Thus the first equation for A'(7r,'i/') takes form 

af (V)(Fo - 27r^) +a±(^)Fi = A'(7r, V)af (V). 

Applying the same procedure for , we get the second equation 

-af (V^)Fr + af (V')(Fo + 2t^^) = A'(7r, V)a±(V'). 
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Finally, we have noticed that the coefficients a^iip) and a^(V') satisfy the system 
of algebraic equations written in the matrix form 

In other words, the main asymptotic corrections to the eigenvalues A.^{r]) are eigen- 
values of the matrix on the left-hand side of (j56p . At the end a simple computation 
provides us the formula 



A'± (tt, ip) = Fo± ^/i^^^^p■^ + \Fi\^ , (57) 

where Fq and Fi are given in (jSOp and (j55p . We emphasise that in the case Fi ^ 
the graphs of the eigenvalues Aj_(7r + ^e^) given by ([3T]) and ([57|) look like the ones 
in fig. 5c. 



3.5. Final formulas 

For the Neumann problem ([1])-([2|) we have Mi = 0, Vi{y) = meas2(w)~2 and 

[/0(x;7r) = e^*^^(meas2(a;))-i 
Substituting these into the formulae ([50|) and (f55|) we obtain 



p 2Q33 + meas3(6l) 2 -^33 + meas3(6') 

i'O = , i-l = TT — . 

meas2(a;j meas2(wj 

By Remark 3.1 it follows that —Q33 > meas3(^) and hence 

Fo<0, |Fi|>0. (58) 

However, for the mixed boundary value problem (j26p 

C/^(x;7r) = e±-^yi(y), 

where {Mi, Vi} is the first eigenpair of the Dirichlet problem (I30p . In this case the 
coefficients Fq and Fi are 

Fo = VyViiyyQ'VyViiy'') + (Mi + tt^Q^s + meas3(^))|Fi(y°)|2 (59) 
Fi = 6-2--" [V,yi(y°)TQV,yi(y°) + (Mi - ^^(^33 _ meas3(^))|yi(y°)n 

- z2vre-2--{Q3if^(y°) +g32^(2/°)}Fi(2/°) 
where Q' is the 2x2 upper left block of the virtual mass matrix Q for the set 

4. Justification of the asymptotics and existence of a spectral gap 

In this section we will first estimate the difference A^(?7) — A^(?7) outside the Ce^- 
neighbourhood of the point r] = it and show that the perturbed eigenvalues A^ (r/) 
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are close to the eigenvalues of the unperturbed problem. This is a rather standard 
step in the asymptotic analysis of eigenvalues in singularly perturbed domains. The 
second step is more elaborated. 

In order to prove the opening of the spectral gap we have to derive approriate 
estimates of the remainder in the asymptotic expansion of A^(7/) constructed in 

7 

section 3 inside the Ce 2 -neighbourhood of the point r/ = vr. To this erid, apply the 



following classical lemma on "near eigenvalues and eigenvectors" (see [271] and also, 
e.g., [H, Ch. 5]). 

Lemma 4.1: Let B be a selfadjoint, positive, and compact operator in Hilhert 
space % with the inner product (•, •)-^. If there exists a number h > and an element 
Wen such that \\W\\h = 1 and \\BW - bW;'n\\ = t £ (0,6), then the segment 
[b — T,b + t] contains at least one eigenvalue ofB. 



4.1. T] far from n 

In this section we will prove the following theorem 

Theorem 4.2 : For any a > there exist C, C{a) > and > 0, such that for 
all e £ [0, ea) 

A°(r?) - C«63- < AUr^) < Al{r^) + Ce^ k = 1,2. 

The following lemma is a direct consequence of the one-dimensional Hardy in- 
equality: if u{b) = 0, then 

b r-b 

\u{r)\'^ dr < i \dru{r)\'^r'^ dr . 
J a 

Lemma 4.3: There exists a positive constant C = C{9,w,x^) such that the fol- 
lowing inequality holds: 

\\{r + e)-\; L\w,)\f < C ^u-L^{w,)f + \\u- {uj \U)f') , 

where r = |x — xo| andU C w is a fixed neighbourhood of the point x^, independent 
ofe. 



We also need an a priori estimate for the eigenfunctions of the problem (j23p - (j25p 
in the singularly perturbed domain zu^. The method developed in [l4] gives such 
estimates in weighted Kondratiev norms 

u';V^{m,) ^ = I ^^^^"'^^^ ^lu\x) ^ dx. (60) 
fc=o -^^^ 

and the step-weighted norms 

\\U^-Vl,{w,)\\ = (||V.f/-y^i(n7,)f + ||C/-y^o„l(z^,)f)^/^ (61) 

For the Neumann problem, the latter norm provides asymptotically sharp estimates 
and we formulate this estimate with references to the general result in P, Ch. 4 
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and 5] and paper [sl] where the detailed explanation of the method is given for a 
concrete boundary value problem. 



A*^. For any /3 G ( — ^, |) the inequality 



Lemma 4.4: Let he an eigenfunction of problem (123p -(l25 p with the eigenvalue 
1 1^ 

'2' 2' 



||[/^y|o(^e)|| <CpK\\U';L\w,)\\ 

is valid while a constant Cp is independent of and e G (0, eg], eq > 0. At the 
same time Cj^ — )■ +oo as (5 ^ 

Proof of Theorem 4.2: 

According to the max-min principle (cf. P, Thm. 10.2.2]) 

Afc ?? = sup mf — — -J- — , 62 

where supremum is taken over all subspaces Ek C H^{we) with co-dimension k — l. 
Every such contains the linear combination of eigenfunctions v = aiUi+a2U2 + 
. . . + akUj^ with |aip + • + la^p = 1, so in the case c^e^ < 1 we conclude 

\\V,v;LH^,W . \\V,v;L\w)f 

^^k\n) S ,, .. r2^„ MI2 - 



\v;L^iw,W ~ \\v;L^rnW-\\v;L^e,W 
^ |aipAO(r?) + |a2pA^(r?) + . . . + |a,pA°(r?) ^ AO(r?) 
~ |ai|2 + |a2p + . . . + lofep - CkC^ ~ I - Cke^ 

<Al{rj){l + Cke^). 

Next we give a proof of the first inequality of the statement for A\. By the 
minimum principle [l|, Thm. 10.2.1], we have 



\\v;L'^{zu) 

Replacing v by the eigenfunction u\ normalised to 1 in L'^{zUe) is not possible 
and we construct an extension of function u\ as follows. Consider the function 
w{e~^{x — x^)) = u\{x) which can be written as 

where is orthogonal to 1 in L'^{d9). Since dO is smooth, function w-^ can be 
extended to a function w-^ G by some fixed continuous linear extension 

operator such that the inequality 

\\V^w^;L\e)\\<C\\w^-H\BR\e)\\ (63) 
<C\\V^w^-L\BR\e)\\=C\\V^w,L\BR\9)\\ 

is valid. Setting u\{x) := ■w-'-{e~^{x — x^)) + wq we derive from ([63|) 

\\V^ul;L\e,)\\ < C\\V,ul;L\BR,\e,)\\. 



October 28, 2011 



2;14 Applicable Analysis BANARU gAPA 



20 

Now we have 

AO. ^^ ||V,.^/f;L2(tu,)f + ||V,.5i;L2(0,)||' 

AVfn) < ^^^5 ^5— ^ 

< Al(??) + ||V,.51;L2(0,)f 
Owing to Lemma 4.4 and the definitions ()60p and ()6ip we have 

where we can take /3 = — ^ + f G ^) ^iid obtain the first estimate. The second 
inequahty for the eigenvalues A^(r/) is proved in a similar manner but with an 
application of the max-min principle. 



4.2. T7 « TT 

Theorem 4.2 ensures the following statement: For any small positive ?7o and ttiq 
one can find eo = e(%, > such that for all e G [0, eo) and 77 G (tt — ryo; + %) 
the interval (Mi + (vr — r/o)^ — mo, Mi + (vr — 779)^ + mo) contains either two simple 
eigenvalues, or one eigenvalue with multiplicity two. 

The immediate objective becomes to justify the asymptotics of the eigenvalues 
constructed in Section 3. 

In Hilbert space 'H'^i'i]) := H}^{'uj^) with scalar product 

we define the operator B^{rj) by the formula 

(S'(?7)u,i;)^.(^) = {u,v)^^, u,v e n^{ri) . (64) 

The operator is self-adjoint, positive and compact. Thus the spectrum of B'^lr]) 
consists of the point /3o = implying the essential spectrum and of the positive 
decreasing sequence of eigenvalues 

/3Kr?)>/3^(r/)>...>/3j(7?)>... (65) 

The variational formulation of the spectral problem (j23p 



gives the following relationship between the eigenvalues of the operator B''{rj) and 
the boundary value problem (j23|) - (p5|) 

= P){ri)-' ■ (66) 
The relationship ([66]) turns the sequence ([65|) into the sequence (fT9]) . 
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We now apply Lemma 4.1 to the operator ()64p . To this end we set 

b± = ii±r' = (Ai(^) + 63A'±(vr,V))-\ (67) 
V±(x) = \\U±;n'{7r + ^e')\\-'U±ix), 

where 

U±{x) = X,{x)\Jlix;7r,iJj) + {1 - X,{x))pI{x^;7t,^P) (68) 
+ Xe{x)W±{x; vr, i/j) + e^X,{x)\J'^{x; vr, ip) + e^U"(x; V), 



and 



Pi(x; ^, V) = Vl (x°; vr, ^) + J] ^(:eO; vr, ^) 

j=i 



1 ^ a^ui, 

+ o a;jXfc- {x'';7r,'ip), 

j,fe=i 

W±(x; TT, V) = eWl{e-\x - x°); vr, V) + e2w|(e-^(x - x°); vr, V) • 

Here the cut-off function is defined by 

X,(x) = l-xe(e-'(x-xO)) 

where xe £ C^([R^) is such that X6» = 1 in the neighbourhood of the set 6. The 
term U" G H^{w) in the sum is added to satisfy the quasi-periodicity condition 
(1251) with 77 = vr^ + e^tp. Due to the definition of the functions Wj., W^, X^, U" 
and xe the function V± belongs to space 7i'^{TT + ipe^). 

To use Lemma 4.1 we have to verify that b± = {£±)~^ and ||[/^; ?^'^(7r + ^/^e^)|| 
are separated from zero, and then we have to estimate 

T± = \\BW± - bW±;'H' {tt + 1^6^)11 

for e > small enough. We proceed with the following assertion 
Lemma 4.5: For all e > there holds the inequality 

\U±;'H'{TT + i;e^)\\-{Mi+TT^y/^\si^y < ce^/^^l + e^/^^;). (69) 

Proof: Since U^(x;7r,^/;) = af{'il>)U^{x) + af{'ip)U^{x), |a='=|c2 = 1, is a solution 
of the problem (f9|)- (lll|) for r] = tt + ipe^, we have 

||V^U^;L2(ro)f = {Mi + tt^). 



The expression on left-hand side of the inequality (I69p is estimated by the sum of 
the expressions 

\\V^lJl;L\e,)\\, 

1 1 ( 1 - X,) (Ui - Pi ; (vr + Ve') 1 1 , 1 1 xe W± ; (vr + ^e^ ) 1 1 , 
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To process the first two norms it is sufficient to observe that the supports of the 
functions V^^e and 1 — axe contained in a ball of radius 0(e), where the Taylor 
formula 



|U^(a;;7r,V) -P±(x;7r,V')| < c|x-a;°|^ < Ce^ (70) 



is valid. 

For the estimation of the norm 1 1 W-t- ; (tt + V-'€'^)|| "wc have to take into ac- 
count the behaviour of the function W± at infinity. A direct calculations leads to 
the estimate 



<ce3(l + e3/2|y,|)2_ 



Finally, using that U'_|_,U" G H^{w), we get the bound for the last two norms 
which proves the statement. □ 



Under the restriction e^/^l'i/'l G (0, V'o] we conclude inequalities 
i±>ce>0, ||W±;H'(7r + ^e=^)|| >Q^ >0, 



with some positive constants q and cu, depending on ipQ, w and 9 only. 
For the estimation of r we use the formula for the norm 



= max I [B^n + ^e3)V± - £±V±, ^)^.(,+^,3) 
<cJ%'max\T±{Z)\. 



Here the maximum is calculated over all functions Z e'H'^{7r + ipe^) with the unit 
norm, and T± stands for 



r±(z) = (v,.ZY±,v,.z)^.(^^)3 

-(A^ + e3A'±(7r,^))(Z^±,Z)i.(^^). 
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T±{Z) can be now represented as a sum of the expressions 

Si\z) = (-(A,. + AO )U0 = 0, 

5i')(Z) = (-[A,,X,](uO -P0),ZU, 

sf{Z) = e(-A,.Wi,xe^)^, - (e5„Wi - V,Ui(xO;7r,V)^n,Z)ae. (71) 

1 ^ S^Ijo 

j,k=i ^ 
sf{Z) = (-[A,,xe](eWi +e'w2 

- (Ai + e3A'±)(eWi + e^W^, xe^)^. 

= e=^(-(A, + A" + A'±)U'± + ^i, 
5i')(Z) = e3(_[A,,X,]U'±,ZU -e6A'±(U'±,X,Z)^^ +e6(v,U",V,Z)^^ 

-e6(A0 +e3^'^)(U",ZU 

Here denotes the last term in equation (j52p . The expressions (j7ip can be 
divided into three types: 

• the terms which appears due to the multiphcation the function U% by the cut-off 
function, 

• the terms which stem from the asymptotic boundary layers and 

• finally, from terms which are asymptotic smooth corrections of U'^. 

To estimate the first type of terms it is enough to observe that the supports of 
V and 1 — Xe are included in a ball of radius 0(e), and to apply the Taylor 
formula (I70p in this ball. Hence the Holder inequality gives us the upper bound 

ce"^\\Z-H\wM 

for first type terms. 

Components of the third type, except for the surface integrals, are evaluated 
in the same manner owing to the additional multiplier e^. The surface integral 
(9nU^,Z)gg can be estimated by the trace inequality and the asymptotic be- 
haviour of the function Uj_ in the vicinity of the point which provides us the 
point wise estimate |U'_|_(x)| < clog |x — 

Now it is suffices to estimate the contribution coming from the asymptotic bound- 
ary layer terms. For this assessment we use the main parts of the functions 
and W^. Let us denote by 

eWi(e; vr, ^) + e^W^^; vr, ^) = 6^° (x) + w±(x). 

The asymptotic behaviour of w*^ in the neighbourhood of x*^ is |x — x*^!"^ and 
w±(x) = O(e^). Thus by Lemma 4.3, or more precisely by the Hardy-type inequal- 
ity, 
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I Ai (vr) (xw± , Z)^^ + e'A'^ (vr, (xw^ , Z)„ J 
< ce^/2|||x-x°|-^^;^2(t:J7,)|| 
<ce'^/^\\Z;H\uj,)\\. 

The evaluation of the remaining terms shall be made at the expense of an addi- 
tional factor e^, or due to the replacement of by w-t leads to an improvement 
in the estimates. 

Collecting obtained inequalities we obtain that r < ce^/^. Thus by Lemma 4.4 
the operator B'^^ir + t/^e^) has just two eigenvalues /3-i-(7r + tpe^) in the interval 
[^ii-ce3/2,£^i + ce3/2]. 

Now the relationship (|66]) of the spectral parameters of the operator B gives the 
result. 

Theorem 4.6 : There exist C > 0, eo > and ipo > such that for all e € (0, eo] 
under the condition \ip\e^^'^ < ■00 the problem (^-([^ with rj = ijj + -ne^ has two 
eigenvalues satisfying the inequality 

|A^(7r + V'e3) - Ml -vr^ -e3A'±(7r,V')| < Ce'^''^ , 



where is taken from (51). 



4.3. The detection of the gap. 

To detect a gap in the spectrum we apply the results of previous computations. 
First, if jry — 7r| > il^oe^/'^ with ■0o > 0, then by means of Theorem 4.2 we will show 
that the interval 

(Ml + vr^ + (Fo - |Fi|)e3 + Mi + vr^ + (Fq + \Fi\)e^ - ce^/^) (72) 

is free of spectrum. Otherwise, in the case — 7r| < V'oc'^^^ we make use of much 
more elaborated asymptotic result of Theorem 4.6 to show that in the interval (I72p 
there is no eigenvalues. This provides us the existence of the spectral gap. 
To realise this scheme, we first observe that, if with some ■0o > 

r? e (0, vr - Voe^^^] U [vr + Voe^''^ 27r) , 

then by Theorem 4.2 

A^r/) < A?(7?) + Ce^ < Ml + TT^ - 2vr|V'oe3/2| + V-ge^ + Ce^ (73) 
and choosing 5 G (0, ^] in Theorem 4.2 

^\{^) > A2(f?) - Ce^-'^ > Ml + vr^ + 2vr|0oe^/^| + ^o^^ - C'e^"" • (74) 

Because of the terms ib2vr^/;o£^ on the right hand side of ()73p and (|74p we may 
choose e > so small that 



Af (??) < Ml + ^2 ^ (Fo - |Fi|)e3 + ce 



7/2 
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and 

K\{r,) > Ml + vr^ + (Fq + \Fi\)€' - ce'"^ , 

and therefore the interval ()72p does not include any eigenvalues. 
If ?7 = TT + ipe^ G [tt — V'oe^^^) + 4'o(-^^'^], then by Theorem 4.6 we have 



and 



< Ml + vr^ + (Fo - \Fi\)e^ + ce^/^ 

A|(7?) > Ml + vr^ + (Fo + vWn;^) - ^^g^ 
> Ml + TT^ + (Fo + |Fi|)e3 - 



so interval (j72p also is free of eigenvalues and we can formulate the main result of 
the paper. 

Theorem 4.7: Assume that \Fi\ > 0. Then there exist positive numbers eg and 
Co such that, for e G (0,eo], the essential spectrum of the problem dl])-© or 
has a gap of length £{e), 

\i{e)-2\Fi\e^\ < 2coe^/2 
located just after the first segment Ti(e). Moreover, the gap includes the interval 
(Ml +Tr^ + (Fo - |Fi|)e3 - coe^/^^Mi + vr^ + (Fq + jFiDe^ + coe^/^) _ 



4.4. Asymptotic formulae for the position and size of the spectral gap 

We proceed with discussing the Neumann problem ([l|)-([2]). Based on the relations 
(j58p we first conclude that the condition |Fi| > of Theorem 4.7 is always fulfilled, 
so that the gap /(e) = (^_(e), A_|_(e)) opens in the spectrum. Moreover, its length 
is 

^^,^3 meas3(g)-Q33 ^Q(^i)^ 
meas2(i^) 

where the endpoints admit the asymptotic form 

^_(e) = + 2vrV^^ + 0(6i), 
meas2(w) 

^^_(6) = .2 + 2vrV^^^^^ + 0(ei). 

meas2((^) 

Since Qas < (see Remark 3.1), the gap covers the point vr^. 

For the mixed boundary value problem ()26p let be the maximum point of 
the positive eigenfunction Vi. Then VyVi{y^) = and the formula (j59p leads to a 
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similar simple situation as for the Neumann problem, because then the number 

e^-"'Fi = V„Fi(9")T(?'V,Vi(/) + (Ml + ,r2(nieass(«) - e33))|Vi(9»)l' 

_2.i(03,F,(,»)^+Q3.V.(!/»)^) 
^ dyi dy2 ' 

= (Ml + 7r2(meas3(^) - Q33)|T4(y°)|' 
is for surely positive. The interval /(e) has the endpoints 

A_{e) = Ml + TT^ + 27rh^Qs^\V,{y^)f + 0(ei), 

A+{€) = Ml + TT^ + 2e3(Mi + 7rVeas3(^))|Fi(/)P + 0(e5) 

and its length is 

2e3(Mi +7r2(meas3(0) - Q33))|V^i(y°)|' + 0(e')- 

The point Mi + tt^ is contained in the interval /(e). 

Theorem 4.7 ensures the opening of the in the case > 0. However, in the 
mixed boundary value prol)lcm varying the position of the point in uj wc may 
come across the relation Fi = 0. In this case the asymptotic formulae obtained 
here do not allow us to conclude the existence of the gap in the spectrum. 

To show that there exist points for which \Fi \ = 0, we assume that Q12 = Q13 = 
0, for example due to the central symmetry of the set 6. Then the expression 

e2--o^^ = V,Fi(y°)Tg'V,Fi(yO) + (Mi + vr2(meas3(^) - g33))|^i(y°)P 

is real. Moreover, if y^ is the maximum point of Vi, then the expression is positive; 
but it becomes negative when y^ locates close to the boundary dui where the 
homogeneous Dirichlet condition is imposed. This is due to the fact that the second 
term, which is positive, diminishes, while the first term remains strictly negative all 
the time according to Remark 3.1, because the normal derivative dnVi is negative 
everywhere on du. Hence there must be a subset in u> where Fi = 0. To prove or 
disprove the existence of a spectral gap one needs to construct higher order terms 
in the asymptotic ansatze. 
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